Conserved charges of order-parameter textures in Dirac systems 
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A simple expression for the induced fermion current in the presence of a texture in mass-order- 
parameters in two-dimensional condensed-matter Dirac systems is derived using the representation 
theory of Clifford algebras. In particular, it is shown that every texture in three mutually anticom- 
muting order parameters, in graphene for example, implies an induced density of a properly defined 
conserved charge. The sufficient condition for the general charge to be the familiar electrical charge 
is that the remaining two anticommuting order parameters allowed by the particle-hole symmetry 
are the two phase components of some superconducting order. This allows eight different types of 
electrically charged textures in graphene or in the 7r-flux Hamiltonian on the square lattice. Gener- 
alized charge of mass-textures on the surfaces of thin films of topological insulators, or in spinless 
Dirac fermions hopping on the honeycomb lattice is also discussed. 



I. INTRODUCTION 

It is well known that in Fermi systems described by 
the minimal Dirac Hamiltonian in two spatial dimensions 
with a general and inhomogeneous mass term, 

H (m) = 170(^171 +P272 + imi + m 2 73 + m 3 7s), (1) 

where m, = mj(af), the induced fermionic current is pro- 
ortional to the conserved topological mass-current: 0- 
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Am • (d u rh x d\fh). 



(2) 



Here the Hermitian four-dimensional 7-matrices satisfy 
the anticommutation relations {7 M ,7„} = 2<^„, and \& = 
^jo- This, for example, implies that a Skyrmion config- 
uration in the mass m in total carries a single unit of the 
Z7(l) gauge charge. Indeed, in graphene, as a condensed 
matter embodiment of the Dirac equation, the Skyrmion 
in the vector order parameter that describes the anoma- 
lous spin Hall state can be understood as two identical 
copies of the above Dirac Hamiltonian, yielding the total 
electrical charge of two. [3, Q Bilayer graphene, the ab- 
sence of Lorentz invariance of its Hamiltonian notwith- 
standing, can also be shown to support Skyrmions in the 
anomalous spin Hall state, but of the twice larger charge 
of four.JBI, |(| This purely topological mechanism of elec- 
tron pairing opens up new and unconventional possibili- 
ties for the emergence of superconductivity in such Dirac 
systems @, HI , and has also been discussed in the context 
of deconfined criticality. [1, H, ^ 

In the single and double layer graphene, however, the 
low-energy Dirac (or more precisely, Dirac-like, for a bi- 
layer) Hamiltonian is eight-dimensional, and if Cooper 
pairing is to be accommodated, sixteen-dimensional ma- 
trix operator. It therefore may not come as a surprise 
that when these Hamiltonians can be block-diagonalized, 
with the elementary blocks of the form in Eq. (1), the 
electrical charges of the individual blocks may mutually 
cancel. This is precisely what happens if the Skyrmion 



is formed out of the three components of the Neel anti- 
ferromagnetic order parameter in the single-layer, Q or 
out of the layer-antifcrromagnetic order parameter in the 
bilayer graphene, § when the residual chiral symmetry 
of the Hamiltonian guarantees vanishing of the electri- 
cal charge. In this sense the anomalous spin Hall state 
may appear as if enjoying a privileged position among 
many possible order parameters in the Dirac electronic 
systems, and being possibly the unique "dual" to super- 
conductivity. [13] 

This raises the obvious question addressed in this pa- 
per: what exactly is the condition on the three mass 
terms (order parameters) being wound, so that the 
Skyrmion, or a more general texture, develops the elec- 
trical charge? The answer follows from a more general 
result, which we discuss first, and prove later in the 
text: if the low-energy excitations in presence of a tex- 
ture in three mutually anticommuting order parameters 
mi, i = 1, 2, 3 are described by the general Dirac- Nambu 
particle-hole-symmetric Hamiltonian 



H = iTaipiTi + p 2 T 2 + imi + m 2 T 3 + m 3 T 5 ) 



(3) 



where the Hermitian T-matrices are d- dimensional (d = 
8,16), and {r^r^} = 25^ v , then the induced fermion 
current is given by the following expression: 



(*r M X*) = — e^ x m ■ (d u m x d\m), (4) 

647T 

with ^ = "I^ To, and with the matrix X as 



x — roTi^rars. 



(5) 



This means that every texture in three mutually anticom- 
muting order parameters, insulating or superconduct- 
ing, carries a distribution of some generalized conserved 
"charge" (^X^), determined by the matrix X, which 
also commutes with the Hamiltonian H. The particle- 
hole symmetry prohibits X to simply be the unit matrix, 
but the general charge reduces to the electrical charge if 
the matrix X happens to be the operator for the particle 
number in the given representation. 
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In graphene, as a representative of a class of two- 
dimensional condensed matter Dirac systems which obey 
parity and the time reversal symmetries for each spin 
component separately, [ll| there are two sublattices, two 
Dirac points (or valleys), and two projections of the elec- 
tron spin, so including the particle-hole doubling, it fol- 
lows that the relevant dimension of the Dirac-Nambu 
Hamiltonian is d = 16. Quite generally for this class of 
systems, the particle-hole symmetry allows 36 different 
mass terms, which can be grouped into 56 fiv etup lets of 
mutually anticommuting order parameters . fl2l [l3l] Let us 
consider the particular, and as it turns out the unique, 
(Appendix C) fivctuplet to which the matrices 1^, i = 

0, 3, 5 belong. One can always choose the particularly 
convenient "Majorana" representation in which the ma- 
trices Ti and T2 are real, and i = 0,3,5,6,7 are 
imaginary, where Tq and T-? are the further two matrices 
that complement To,r3,r5 to the aforementioned five- 
tuplet. [12], |l4j In this representation the Hamiltonian 
H is imaginary, and the (antilinear) particle-hole sym- 
metry, under which H changes sign, is just the complex 
conjugation. It will then follow that 

x = tvf 7 , (6) 

1. e. the matrix X is the generator of "rotations" be- 
tween the remaining two matrices in the fivetuplet, F@ 
and IV If X is the number operator, the latter two must 
stand for the two phase components of a superconduct- 
ing order parameter. Since in graphene- like systems there 
are four possible types of the superconducting mass-gaps, 

15] and each pair of anticommuting mass-matrices be- 
longs to two different fivetuplets, [14[ (discounting the 
rotations of electron spin) there are precisely eight dif- 
ferent types of triplets of order parameters which can be 
combined into an electrically charged texture. We discuss 
some examples of these in the concluding section. 

Although the matrix X, being imaginary in Majorana 
representation, respects the particle-hole symmetry, it 
is not the unique such matrix that commutes with the 
Hamiltonian H. There are in fact three more, as dic- 
tated by the quaternionic nature of the real representa- 
tion of the Clifford algebra C(2, 5) [l6[ furnished by the 
set of sixteen-dimensional matrices {Ti,iTj}, i = 1,2, 
j = 0, 3, 5, 6, 7. [3] Therefore, the conservation of the 
particle-number in the Dirac-Nambu Hamiltonian alone 
is a necessary but not also a sufficient condition for the 
finite electrical charge of the texture. The majority of the 
insulating triplets in graphene will in fact yield precisely 
zero density of the electrical charge. 

It is also interesting to consider the case when the di- 
mension of the matrices in the Hamiltonian H is d = 8, 
describing the spinless fermions hopping on honeycomb 
lattice, or the two boundary surfaces of a thin film of 
topological insulators, once Nambu's particle-hole dou- 
bling is implemented. In this case there are only 10 
possible mass terms, and 6 different triplets of mutu- 
ally anticommuting masses. (Appendix A) Since again 
Ti are real for i = 1,2 and imaginary for i — 0,3,5, in 



this case there are no further matrices that would anti- 
commute with the entire set of To,^ i — 1,2,3,5. [l4| 
The matrix X is then the unique matrix which respects 
the particle-hole symmetry (that is, which is imaginary 
in the Majorana representation) and that also commutes 
with the Dirac-Nambu Hamiltonian. In this case, and in 
contrast to the case when d = 16, it is not only necessary 
but also sufficient that the Hamiltonian H conserves the 
particle number, in order for the texture to be electrically 
charged. In fact, the electrically charged triplet of order 
parameters is in this case unique. 

In the following two sections we prove the Eqs. (4) 
and (5) for the matrix dimensions of d = 16 and d = 8, 
which involve different types of representations of Clif- 
ford algebras and require related, but somewhat differ- 
ent proofs. A reader not interested in the mathematical 
details may skip over these two sections and go directly 
to the Section IV. In Appendices A and B we provide 
some simple (and for our purposes helpful) facts about 
the requisite Clifford algebras and their relevant represen- 
tations. In Appendix C we prove our starting assertion 
that any triplet of anticommuting masses in the sixteen- 
dimensional Dirac-Nambu Hamiltonian is always a subset 
of some fivetuplet of compatible (anticommuting) order 
parameters. 

II. GRAPHENE-LIKE SYSTEMS: d = 16 

Let us now prove the Eq. (4). We start with 
the case relevant to the electrons in graphene, or, 
more generally, to any other Dirac Hamiltonian deriving 
from the spin-independent electron hopping on a two- 
dimensional lattice and respecting parity and time rever- 
sal, when the dimension of the matrices is d = 16.[li1| 
Affleck-Marston 7r-flux Hamiltonian on the square lat- 
tice, for example, falls into the same class. The ma- 
trices Ti, F2, iFo, iF3, iTg, iT^, iT? are then real in Majo- 
rana representation, and form an irreducible real repre- 
sentation of the Clifford algebra C(2,5). 14)] The orig- 
inal matrices T\, F2, Fo, F3, F5, F6, F7, however, may al- 
ternatively be viewed as forming a particular sixteen- 
dimensional complex representation of C(7,0). Let us 
denote this representation by Repie(2,5), as a reminder 
that the representation is 16-dimensional, and that two 
matrices in it are real, and the remaining five are imag- 
inary. On the other hand, it is well known that any 
Clifford algebra C(2n + 1,0) has only two inequivalent 
irreducible complex representations which are both 2"- 
dimensional, and which differ only in the sign of an odd 
number of matrices. [17] For C(7,0), these will be, m our 
notation, Irep%(A : 3) and 7repg(4, 3). It then follows that 
either 

Rep 16 (2, 5) = Irep 8 (4, 3) © Irep 8 (i, 3), (7) 

or 

Rep 16 {2, 5) = Irep' s {4, 3) © Irep' s (A, 3). (8) 
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In other words, the matrices ri, T2, Tq, r 3 , Tq, 1^ can 
be block-diagonalized into an orthogonal sum of two 
equivalent irreducible eight-dimensional representations 
of C(7, 0). This is in accord with the fact that C(2,5) 
has two inequivalent real 16-dimensional representations, 
which correspond to the two possibilities in Eqs. (7) and 
(8). This result follows immediately from the observation 
that the only remaining possibility, namely that 

Rep 16 {2, 5) = Ire Ps (A, 3) © Irep' s (4, 3), (9) 

i.e. that r p = ct 3 C8 a M , where {a^,a^} — 25^ v and 
a-matrices form an eight-dimensional irreducible repre- 
sentation of (7(7,0), would imply that there are only 
two matrices which commute with RepiQ(2,5): <jq ® Is 
and 03 ® 1%, where Is is the eight-dimensional unit ma- 
trix. This is because the unit matrix Is is, by Schur's 
lemma, the only matrix that commutes with all seven of 
a-matrices. This would contradict, however, the fact that 
the sixteen-dimensional real representation of C(2,5) is 
quaternionic, and that consequently RepiQ(2,5) actually 
has four Casimir operators. [14], Qjl HH It is easy to see 
that both the decompositions in Eqs. (7) and (8) con- 
form to this fact, and these are simply, after the block- 
diagonalization, cto (S> Is, and a <8 Is- 

By a judicious unitary transformation the Dirac- 
Nambu Hamiltonian H can therefore be transformed into 
an orthogonal sum of two identical eight-dimensional 
copies: 

H — CT ® iao(p\a\ + p2a 2 + im\ + m 2 a 3 + m 3 a 5 ), (10) 
with the five eight-dimensional a-matrices representing 
I 



with the two factors 01 ® 17172 and cti (8)17172 as precisely 
the two remaining matrices that anticommute with the 
set 

CT <8«7o7i> CT o®«7o72,cr 3 ®7o,CT3<g)i7o73,CT3<g)i7o75, (15) 

which appears in the second factor in Eq. (11) when 
written explicitly. Reversing the unitary transforma- 
tions used to bring the Hamiltonian into the final block- 
diagonal form in Eq. (11), Eq. (14) implies Eq. (6), or 
equivalently, Eq. (5). 

III. TOPOLOGICAL INSULATOR THIN FILM: 

d = 8 

One can similarly establish the validity of Eq. (4) 
when the dimension of the matrices in the Hamiltonian 
is d = 8, with the details of the proof being only slightly 



C(5,0). But C(5, 0) also has only two inequivalent irre- 
ducible representations, which are then four dimensional. 
Again, the very existence of two further anticommuting 
matrices T§ and I\ guarantees that the particular matri- 
ces «0j Qij cti, 0L3, 05 appearing in the above decomposi- 
tion must now be an orthogonal sum of two inequivalent 
representations of C(5, 0). A more detailed proof of the 
last statement will be given in the Section III when we 
discuss the case of d = 8. 

In sum, any sixteen-dimensional Dirac-Nambu Hamil- 
tonian H, obeying particle-hole symmetry by construc- 
tion, can be transformed into 

H = ct © (H (m) © H (-m)) (11) 

where Hq is the irreducible four-dimensional Hamiltonian 
in Eq. (1). 

It then becomes evident that 

(#r M (<7 ®<To®i4)tf> =0 (12) 

always, because of the cancelation of the two contribu- 
tions with opposite signs. On the other hand, in the 
block-diagonal representation in Eq. (11), 

(*r M (£r ®<r 3 = —epvxm-idvmxdxm), (13) 
since the sum of the two contributions now add up. The 
factor d/8 = 2 accounts for the doubling due to spin-1/2 
degree of freedom, but not for the particle-hole doubling. 
The inclusion of the latter would, of course, lead to over- 
counting. We can also recognize the matrix appearing in 
the last expression as the matrix X in Eq. (4): 



(14) 

I 

different. The matrices I^, T2, iTo, iTa, 1T5 then provide 
a real eight-dimensional representation of the Clifford al- 
gebra C(2,3). (Appendix B) The original five matrices 
ri, T2, To, r§, on the other hand, can also be viewed 
as a particular complex representation of C(5,0); call 
it Reps(2,3). As mentioned earlier, C(5,0) has only 
two irreducible complex representations; let us call them 
Irep4(3, 2) and Irep' 4 (3, 2). One can then show that 

Rep 8 (2, 3) = Irep 4 (3, 2) © Irep' 4 (3, 2), (16) 

i. e. Reps(2,3) can be decomposed into an orthogo- 
nal sum of two inequivalent irreducible representations 
of C(5,0). Again, this follows from realizing that the 
remaining possibility that Reps(2,3) is an orthogonal 
sum of two equivalent representations of C(5,0) would 
imply that the matrices can be transformed into the 
form = cto ® 7^. This, however, would mean that 
there is no matrix that anticommutes with all the five 



X = ctq (8) CT3 <g> I4 = i((JQ <8 ct 2 (8 i7i72)(CTo (8 cti (8 47172), 
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matrices in i?epg(2,3). This would contradict the fact 
that i?ep 8 (2, 3) can be obtained by simply dropping two 
real matrices in 7reps(4, 3), so that two further matrices 
which anticommutc with the entire representation must 
exist. (This construction is discussed in Appendix B.) 
Indeed, from the correct decomposition in Eq. (16) we 
find these to be o\ ® I4 and 02 <E) I a- 

The decomposition in Eq. (16) also means that the 
eight-dimensional real representation of the Clifford al- 
gebra C(2,3), up to a unitary transformation of course, 
is unique. 



i. e. as precisely the product of the five matrices which 
appear in the block-diagonal form of H in Eq. (17), im- 
plies the Eq. (5) in the original representation. 



IV. DISCUSSION 

Although formally the expression in Eq. (4) appears 
the same, the condition for the topological current to be 
the electrical current when d = 8 and when d — 16 is 
quite different. One can understand the principal differ- 
ence in somewhat simpler terms as follows. When d = 8, 
the Dirac-Nambu Hamiltonian is formed by the usual 
Nambu's particle-hole doubling of some original d = 4 
Dirac Hamiltonian which described only particles. The 
mass-terms in this original four-dimensional representa- 
tion stand only for the insulating order parameters, and 
the possible superconducting order parameters are ex- 
cluded. Since the dimension of matrices is then t! = 4we 
can maximally have five mutually anticommuting matri- 
ces, and thus there is a unique triplet of insulating masses 
that can form the textures under consideration. The orig- 
inal result in Eq. (1) then directly applies. For spinless 
fermions on honeycomb lattice, the triplet in question 
consists of two Kekule bond density waves and the charge 
density wave. For two surfaces of a topological insulator 
slab, these are the two components of the excitonic in- 
sulator [l9| and the surface magnetization. One can also 
show that when d = 8 there exist three more supercon- 
ducting mass order parameters, each contributing two 
linearly independent mass-matrices. (Appendix A) Fi- 
nally, there is also the quantum anomalous Hall state, 
bringing the total number of possible mass-order param- 
eters to ten. The latter matrix, however, commutes with 
all of the remaining nine, and thus cannot be used to 
form the textures we are discussing. The product of all 
four-dimensional matrices that appear in Eq. (1) is the 
unit matrix, so Eq. (2) can also be understood as the 
special case of the Eq. (5), modulo the prefactor in the 



When d = 8 we can therefore transform the Dirac- 
Nambu Hamiltonian into 

H = H (m) ®H {-m), (17) 

and therefore the induced fermionic current is 

(*iyo- 3 ® I 4 )#) = —e^m • (d v fh x d x m). (18) 
Again, recognizing the sought matrix X in the above as 



(19) 

I 

latter equation. 

When d = 16, the original pre-particle-hole-doubled 
Dirac Hamiltonian is eight-dimensional, and thus twice 
the size of the irreducible representation of the five mu- 
tually anticommuting matrices, as in Hamiltonian (1). It 
can be then block-diagonalized either into an orthogonal 
sum of either two equivalent, or two inequivalent repre- 
sentation of C(5,0). In the latter case, one can always 
find two matrices which anticommute with the original 
eight-dimensional Dirac Hamiltonian. Since only the in- 
sulating order parameters are discernible in the original 
d — 8 representation, this means that the original triplet 
of insulating orders belongs to the fivetuplet of order pa- 
rameters which are all insulating. The product of the five 
matrices in the Dirac Hamiltonian is then the generalized 
charge of the texture, which differs from the electrical 
charge. In the former case, in contrast, the matrices an- 
ticommuting with the Dirac Hamiltonian do not exist. 
That precisely means that upon particle-hole doubling 
of the d — 8 Dirac Hamiltonian into the d = 16 Dirac- 
Nambu version of itself, two such anticommuting matri- 
ces, which are then guaranteed to exist, would have to be 
superconducting order parameters. The product of the 
five matrices in the Dirac Hamiltonian is then unity, and 
the generalized charge is precisely the electrical charge. 

When d = 16, we can finally view the problem along 
the lines of the ref. [l4[ as well. Assuming two mutually 
anticommuting insulating mass-order-parameters, there 
are two different mutually exclusive ways to complement 
them to the maximal number of five. |14j There are thus 
six remaining choices for the third mass that can be com- 
bined with the first two into a texture. One fivetuplet to 
which the two chosen insulators belong is completely in- 
sulating. Choosing the third order parameter from that 
set leads to the non-electrical conserved charge of the 
texture. The other set has one insulating and two super- 
conducting orders, in addition to the two insulators we 
began with. Choosing the third matrix to be that single 
remaining insulator in the mixed fivetuplet endows the 



I 

X = 0-3 ® I 4 = (cr ® «7o7i)(co (8) £7072) (^3 ® 7oX°3 (g) 17073X03 ® «7o75), 
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texture with the electrical charge. 

Consider, for example, the x- and y-components of the 
Neel order parameter in graphene. |20j If the third order 
parameter in the texture is the remaining z-component 
of the same Neel order parameter, the remaining two 
orders in the fivetuplet are the two Kekule (singlet) 
bond-density waves, O Hl[ and the generalized charge 
is the "valley spin". If the third component in the 
texture is the z-component of the anomalous spin Hall 
state, on the other hand, the remaining two are the two 
(phase) components of the z-component of the p-wave 
superconductor. [TM |22| The generalized conserved charge 
is then just the electrical charge. 

If in the previous example the texture is formed out of 
the two phase components of the z-component of the p- 
wave superconductor and the z-component of the anoma- 
lous spin Hall state, in turn, the generalized charge will 
be proportional to the z-component of the conventional 
electron spin, which rotates the x- and y- components of 
the Neel order parameter into each other. Eq. (4) would 
yield the total value of the generalized charge to be two; 
since the spin is (53) = X^) /2, the value of the total 
integrated third component of the spin correctly comes 
out as unity. 



8. Consider, the Hamiltonian in the Majorana represen- 
tation: 

H = p 1 R 1 +p 2 R 2l (20) 

where R\ and R 2 are two real, mutually anticommuting, 
eight-dimensional Hcrmitian matrices. H is evidently 
imaginary, and the particle-hole transformation is the op- 
eration of complex conjugation. 

We can consider always R\ and R 2 as a part of max- 
imally seven mutually anticommuting eigth-dimensional 
Hcrmitian matrices, which we assume all square to unity: 

{Ri, R2, R3, R4, h, h, h}, (21) 

where Ri are real, and Ii imaginary. This is an irre- 
ducible complex representation of C(7, 0), and i?4 = 
iRiR 2 R 3 I\I 2 I 3 . Given Ri and R 2 , one can then form 
ten linearly independent imaginary Hermitian matrices 
that anticommute with both. We organize them as fol- 
lows: 

Vi = (Ii,h,h), (22) 
V 2 = {iR 3 hh,iR 3 hh,iR 3 hhl (23) 



V. CONCLUSION 

We showed that in graphene-likc two-dimensional 
Dirac systems one can always define a conserved 
fermionic current induced by a texture in the three anti- 
commuting order parameters, and derived the sufficient 
condition for the conserved charge to be the electrical 
charge. Our reasoning was based on the properties of 
the relevant representations of certain Clifford algebras 
closed by the matrices appearing in the particle-hole sym- 
metric Dirac-Nambu Hamiltonian with a mass-texture. 
In particular, from the fact that any three anticommuting 
order parameters can always be embedded into a larger 
fivetuplet of such anticommuting order parameters, we 
derived the expression for the conserved current, Eqs. 
(4) and (5). The sufficient condition for the electrical 
charge of the texture is that the two missing order pa- 
rameters in the fivetuplet are the two phase components 
of some superconducting order. 
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VII. APPENDIX A: TEN MASSES FOR d = 8 

We here discuss the masses of the Dirac-Nambu Hamil- 
tonian allowed by the particle-hole symmetry when d = 



% = {iRihh^RJxh.iRihh), (24) 

and 

S = iR\R 2 . (25) 

The last matrix represents the anomalous quantum Hall 
state which commutes with the others: 

[S,V ij } = 0. (26) 

Among the remaining nine, 

{V m k, V n k] — {Vkm, Vkn} = 2S mn , (27) 

for k — 1,2,3, whereas all other pairs commute. There 
are thus only six possible triplets of mutually anticom- 
muting mass terms that can form the textures under con- 
sideration: 

3 

^mft, (28) 

i=l 

or 

3 

^m,Vt, (29) 

k = 1, 2, 3. Out of these six, one and only one will con- 
sist of exclusively insulating order parameters, and thus 
induce the electrical current. 
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VIII. APPENDIX B: ALMOST COMPLEX 
REPRESENTATION OF (7(2, 3) 

One can also understand some of the features of the 
eight-dimensional real representation of C(2, 3) used in 
the text. Namely, from the seven matrices forming the 
representation of C(7, 0) in Eq. (21), we can form 

{R u R2,iIuiI 2 ,iI 3 } (30) 

as the eight-dimensional real representation of the Clif- 
ford algebra C(2,g) with the maximal value of q, which 
is q = 3. This representation is "almost complex" in the 
terminology of Okubo 17] since besides the unit matrix 
it allows one more real matrix that commutes with all 
the matrices in the representation; it is 

R 3 R 4 . (31) 

One can then also observe that Reps(2, 3), given by Eq. 
(30) without the factors of the imaginary unit "i" in the 
last three terms, necessarily allows for two further an- 
ticommuting matrices, as claimed and used in the text: 
these are i?3 and R4. 

IX. APPENDIX C: EMBEDDING OF A 
TRIPLET INTO A FIVETUPLET FOR d = 16 

The starting point of this work was that given two real 
sixteen-dimensional anticommuting matrices T\ and T2, 
and three anticommuting imaginary matrices Tq, L3, and 



L5, one can always find two further imaginary matrices 
r 6 and Tf, so that any chosen pair out of the seven ma- 
trices anticommutes. 

The proof follows from observing that the set {Ti, iTj}, 
i = 1,2 and j — 0,3,5 represents a sixteen-dimensional 
real reepresentation of the Clifford algebra C(2,3). 
C(2, 3), on the other hand, has a unique real irreducible 
representation which is eight dimensional. fl~7l | We saw 
this already as following from the decomposition in Eq. 
(16), for example. Therefore, by a unitary transforma- 
tion we can write 

r, = ct ® a, (32) 

where a, are eight-dimensional, and real for i = 1, 2, and 
imaginary for i = 0,3,5. Since in the eight-dimensional 
space there exist then two additional linearly indepen- 
dent real matrices, call them a$ and a-j (Appendix B), 
two additional imaginary anticommuting matrices are, 

r 6 =cr 2 ®a 6 , (33) 

and 

r 7 = 02 ® a 7 . (34) 

The pair is obviously unique. Three non-trivial Casimir 
operators for the real representation of C(2, 5) given by 
{Ti, iTj}, i = 1, 2 and j — 0, 3, 5, 6, 7 are then u\ ® a%a 7 , 
03 £3> aga7, and i<J2 ® Is- Together with the unit matrix 
they close the quaternionic algebra. 
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